Geospatial Optimization Problems 
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Abstract — There are numerous applications which require the 
ability to take certain actions (e.g. distribute money, medicines, 
people etc.) over a geographic region. A disaster relief organiza- 
tion must allocate people and supplies to parts of a region after 
a disaster. A public health organization must allocate limited 
vaccine to people across a region. In both cases, the organization 
is trying to optimize something (e.g. minimize expected number 
of people with a disease). We introduce "geospatial optimization 
problems" (GOPs) where an organization has limited resources 
and budget to take actions in a geographic area. The actions result 
in one or more properties changing for one or more locations. 
There are also certain constraints on the combinations of actions 
that can be taken. We study two types of GOPs - goal-based 
and benefit-maximizing (GBGOP and BMGOP respectively). 
A GBGOP ensures that certain properties must be true at 
specified locations after the actions are taken while a BMGOP 
optimizes a linear benefit function. We show both problems to 
be NP-hard (with membership in NP for the associated decision 
problems). Additionally, we prove limits on approximation for 
both problems. We present integer programs for both GOPs that 
provide exact solutions. We also correctly reduce the number of 
variables in for the GBGOP integer constraints. For BMGOP, we 
present the BMGOP-Compute algorithm that runs in PTIME 
and provides a reasonable approximation guarantee in most 
cases. 

I. Introduction 

As geo-located social network data becomes more common 
with sites such as FourSquar^H and programs such as Real- 
ityMinin£0, it becomes desirable to reason about such data. 
There are numerous applications which require the ability to 
take certain actions (e.g. distribute money, medicines, people 
etc.) over a geographic region. For instance, a disaster relief 
organization must allocate people and supplies in a region 
after a disaster. A public health organization needs to allocate 
limited vaccine stocks to people across the region. A govern- 
ment needs to allocate funds for education or unemployment 
training across a region. However, allocating any resource will 
cause certain effects - some desirable, some not - based on 
the network connections among geographic locations. In this 
paper we present a formal framework that allows reasoning 
about such geo-located data in order to answer certain queries 
where we have some desired goal to achieve as the result of 

1 https://foursquare.com/ 

2 http://re alitycommons.media.mit.edu/ 
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J^High-cost area (hi_cost) Group 1 (grp x ) 1 Influential center for group 1 (hqj) 

J^sNon-populated area (non_pop) Group 2 (grp 2 ) 2 Influential center for group 2 (hq 2 ) 




Fig. 1. Locations in a district - contingency groups and unpopulated areas. 



our geographically-based resource allocation - all the while 
considering the complex interactions among locations. 

Figure Q] shows a 2-dimensional map of a region. A political 
candidate can only make so many campaign stops and public 
appeals. We assume that a map Ai is discrete (this is a 
common assumption in most GIS systems) and has coordinates 
drawn from [0, . . . , M] x [0, ... AT] where the bottom left 
corner of the map is the point (0, 0). The candidate wants to 
identify the best places to campaign or make public appeals 
to maximize his exposure. Additionally, the map shows un- 
populated areas, areas where campaigning costs are high, and 
areas dominated by one of two constituent groups. All of these 
factors may affect the set of locations the candidate selects to 
optimize his exposure. 

In this paper, we introduce geographic optimization prob- 
lems or GOPs that capture and solve problems such as 
those mentioned above. This framework allows one to more 
prudently position resources in a manner to achieve a goal 
while considering the complex interactions between locations 
(that may be modeled as a network). The organization and 
contribution of the paper is as follows. Section UT1 formally de- 
fines GOPs - specifically we introduce goal-based and benefit- 
maximizing GOPs (GBGOP and BMGOP respectively). Sec- 
tion |m] shows that both GBGOP and BMGOP are NP-hard 
(with the associated decision problems in the complexity class 



NP). Additionally, we prove non-trivial theoretical limits on 
approximation: if GBGOP were to be approximated within 
the logarithm of the input then NP would have a slightly 
super-polynomial oracle. BMGOP cannot be approximated 
within a guaranteed factor greater than 0.63 unless P=NP. 
Section [IV] presents integer programs to solve both GBGOP 
and BMGOP using an IP solver like CPLEX. In Section [Vj 
we show how to correctly reduce the number of variables 
in the integer constraints for GBGOP. We then develop the 
BMGOP-Compute algorithm in Section [VI] that can quickly 
approximate a BMGOP in polynomial time and provides an 
approximation guarantee. 

II. GOPs Formalized 

Throughout this paper, we assume that KA = [0, . . . , M] x 
[0, . . . , N) is an arbitrary, but fixed "map". We define a logical 
language L whose constant symbols are members of A4 and 
that has an infinite set C var of variable symbols disjoint from 
M. C has a set Q = {g\, . . . , g n } of unary predicate symbols. 
As usual, a term is either a constant symbol or variable symbol. 
If t is a term, then gi(t) is an atom. If t is a constant, then 
gi(t) is ground. Intuitively, if p G A4, then gt(p) says that 
point p has property gi. We use Be to denote the set of all 
ground atoms. Well-formed formulas (wffs) are defined in the 
usual way. (i) Every atom is a wff. (ii) If F, G are wffs, then 
so are F A G, F V G, ->F are all wffs. 

Example 2.1: Consider the map M. cpgn in Figure Q] with 
predicates Q = {hi_cost,non_pop,grp\,grp2,hq\,hq2}. 
The predicate exposure not depicted in the figure cor- 
responds to a candidate receiving exposure in a cer- 
tain area. hi_cost((l, 9)), /igi((4, 3)), non_pop((&, 1)), and 
grp2((5,8)) are all examples of ground atoms. 

A state is any subset of Be- We use S to denote the set of 
all states. Satisfaction of formulas is defined in the obvious 
way. State s satisfies a ground atom A, denoted s \= A, iff 
A G s. s \= F V G iff s \= F or s \= G. s (= F A G iff 
s \= F and s j= G. s \= -*F iff s does not satisfy F. 

Example 2.2: The shading shown in Figure Q] defines 
a state. For example, hi_cost((\, 9)) G s cpg n while 
expos we ((1,9)) £ s cpgn . 

An action maps points to sets of ground atoms. 

Definition 2.1 (Action): An action is a mapping a : A4 — >• 
2 Bc . We use A to denote the set of actions. An action-point 
pair is any member of A x M. 

An action-point pair (a,p) is executed if action a takes place 
at point p. Thus, one can think of (a,p) as saying that action 
a occurs at point p. The result of executing a set SOL of 
action-point pairs in state sq is denoted appl(SOL, sq) and 
is the set (sq U {a(p) \ (a,p) G SOL}). 

Example 2.3: Continuing with example 12.61 our candidate 
has actions A cpgn = {nor, appeali, appeal} where nor 
refers to a normal campaign stop and appeal i, appeal refer 
to public appeals to constituent groups 1 and 2 respectively. 
The actions map to ground atoms as follows. 

nor(p) = {exposure(p')\ -inon_pop(p') A d(p,p') < 1} 
appeak(p) — {exposure(p')\ hqi(p) A grpi(p')} 



The first action says that when a normal campign stop is made 
at point p and p' is a populated place one distance unit or less 
from p, then the candidate has exposure at place p' as well. 
The second action says that if the candidate makes an appeal 
(action) at point p and p is the headquarters of interest group 
grpi, then the candidate has obtained exposure in all places 
associated with interest group grpi. 

Definition 2.2 (Cost Function): A cost function, C : A x 
M -> [0,1]. 

Throughout this paper, we assume the cost function is arbitrary 
but fixed and can be computed in constant time. We also 
assume that if A x AA — {(ai,pi), . . . ,(a m ,p m )}, then Ci 
is used to denote C(ai,pi). 

Example 2.4: The cost function for our example is C^ gn 
and is defined (based on some state s) as follows: 
Ccp gn (a,p) = 1 if hi_cost(p) G s and 0.5 otherwise. 

We also assume the existence of a set of integrity constraints 
IC that specify that certain actions cannot be jointly taken if 
some conditions hold w.r.t. the state — such constraints were 
defined before by 0]. 

Definition 2.3 (Integrity Constraint): If $ is a set of action- 
point pairs and x i s a w ff> tnen ^ ^ X i s an integrity 
constraint. 

When $ ^ x is ground, this says that if x is true, then only 
one action-point pair in $ may be executed. Formally, suppose 
s is a state and is a set of action-point pairs and ^p^is 
ground, (s, $') (= $ x iff either s Y= \ or s \= X an d |3> l~l 
<£>'| < 1. (s, $') satisfies an integrity constraint iff it satisfies 
all ground instances of it. (s, $') |= IC where IC is a set 
of integrity constraints iff (s, $') satisfies every constraint in 
that set. Given a state s and set IC of integrity constraints, we 
use IC S to denote the set of all ground instances of integrity 
constraints in IC where the associated wff x is satisfied by 

Example 2.5: Continuing Example 12.41 let IC cpgn be 
{{appeah({4:,3)),appeal2((10J))} ^ TRUE}. This con- 
straint says that an appeal can be made to either group 1 
or group 2 at their center of influence, but not both — for 
instance, these two groups may have opposing views. 

We now introduce the goal-based geospatial optimization 
problem (GBGOP). This problem takes as input a map M., 
initial state So, set of actions A, cost function C, integrity 
constraints IC, positive real number C, and disjoint sets 
Qi n ,Qout Q Be- Intuitively, C restricts the total cost and 
Qi n (resp. 8 out) is a set of atoms that must be true (resp. 
false) after the actions are applied. Our optimality criteria for 
a GBGOP is to minimize the cardinality of the action-point 
pairs. A GBGOP can be viewed as an abductive inference 
problem (i.e. find a set of actions that lead to the current 
state) - where minimal cardinality is a common parsimony 
requirement. 

Definition 2.4 (GBGOP Solution, Optimal Solution): A 
solution to a GBGOP (M, s , A, C, IC, C, 0;„, Q ut) is a 
set SOL C A x M such that: (i) S( aijPi ) eSOL c i < C, 

formally, IC S = {(* ^ X ) G IC\s \= x} 



(ii) (s , SOL) \= IC, and (iii) appl(s ,SOL) |= 

A solution SOL is optimal iff there is no other solution 
SOL' such that \SOL'\ < \SOL\. 

Our next type of problem is a benefit-maximizing geospatial 
optimization problem (BMGOP) that also considers a benefit 
function, defined as follows. 

Definition 2.5 (Benefit Function): The benefit function, 
B : Be — > 3f£ + maps atoms to positive real numbers. 

Example 2.6: In our running example, we use the benefit 
function B cpgn where B cpgn (A) = 1 if A has the form 
exposureQ and otherwise. 

As with cost, we assume the benefit function to be arbitrary 
but fixed and computable in constant time. We also assume that 
if Be = {Ax, . . . , A n }, then B(A,) is denoted b,. A BMGOP 
takes as input, A4, So, A, C, IC, and C - all defined the same 
as for a GBGOP. Additionally it takes benefit function B and 
natural number k. Here k is a bound on the number of actions 
the agent can take as we attempt to maximize benefit as an 
optimality criteria. 

Definition 2.6 (BMGOP Solution, Optimal Solution): 
A solution to a BMGOP (M, s , B, A, C, IC, k, c) is a 
set SOL C A x M such that: (i) \SOL\ < k and (ii) 
S( 0i> p j ) e soLC i < C, and (iii) (s ,SOL) |= IC. 

A solution SOL is optimal iff there is no other solution 

SOL' SUCh that YsAi€appl(SOL,3o) bl < ^-<Aieappl(SOL',s ) bl - 

III. Complexity Results 

Here, we provide complexity results for GBGOPs and 
BMGOPs. First, we establish both as being at least NP-hard. 

Theorem 1: Given GBGOP (M,s , A, C, IC, C, 6 in , 
Oemt), finding an optimal solution SOL C _4x Ai is NP-hard. 
This result holds even if for each a £ A,p € .M, it is the case 
that yg'(p') € a(p), p' — P - i-e- each action only affects the 
point is is applied to. 

Proof Sketch. We embed the known NP-hard problem of 
SET-COVER [2 1 which takes as input a set of n elements, 
S and a family of m subsets of S, "H = {-Hi, • • • , H m }, 
and outputs %' C. % s.t. the union of the subsets covers 
all elements in S and T-L' is of minimal cardinality. We 
encode this problem into a GBGOP as follows: we set 
Q = {gx, ■ ■ ■ ,g n } - eacn predicate in Q corresponds to 
an element in S, the map, M. consists of a single point, 
p, the actions A — {ax, ...,a m } s..t each action a,iA 
corresponds to an element in H and each is defined as 
follows: di(p) — {J x eH ,{gj(p)}. The cost function C returns 
1 for each action-point pair, 6 in = \j gie g{gi(p)}, <dout = 0, 
and finally, we set sq = 0, IC = 0, C = n. □ 

Theorem 2: Given BMGOP (M,s , B, A, C, IC, k, c), 
finding an optimal solution SOL C A is NP-hard. This result 
holds even if for each a € A,p € M, it is the case that 
Vg'(p') G a(p), p' = p - i.e. each action only affects the point 
is is applied to). 

Proof Sketch. We embed the known NP-hard problem of 
MAX-K-COVER |2] which takes as input a set of n elements, 



S and a family of m subsets of S, H = {Hi, . . . , H m }, and 
positive integer K and outputs < K subsets from T-L s.t. the 
union of the subsets covers a maximal number of elements 
in S. We encode this problem into a BMGOP as follows: 
we set Q = {gx, . . . ,g n } - each predicate in Q corresponds 
to an element in S, the map, Ai consists of a single point, 
p, the function B returns 1 for each ground atom, the set 
A = {ax, ■ ■ ■ ,a m } is set s.t. each action in A corresponds 
to an element in H and each at is defined as follows: 
a i(p) = [J x eH {dj(p)}- The cost function C returns 1 for 
each action-point pair, and finally, we set sq = 0, IC = 0, 
k = K,C = K. □ 

One may think that one can solve GOPs efficiently in prac- 
tice by using fully polynomial time approximation schemes 
(FPTAS). However, by the nature of our constructions used in 
the NP-hardness results, this is not possible for either type of 
GOP under accepted theoretical assumptions. 

Theorem 3: If for some e > 0, there is a PTIME algorithm 
to approximate GBGOP within (1 - e) • ln(|^4 x M\), then 
NP c TIME(\A x M\° {lglglAxMn ) (NP has a slightly 
super-polynomial algorithm). 

Follows from Theorem Q] and [2. Theorem 4.4]. □ 

Theorem 4: Finding an optimal solution to BMGOP cannot 
be approximated in PTIME within a ratio of + e (approx. 
0.63) for some e > (where e is the inverse of the natural 
log) unless P=NP, even when IC = 0. 

Follows from Theorem [2] and [2. Theorem 5.3]. □ 

Next, under some reasonable assumptions, the decision 
problems for GBGOP/BMGOP are in-NP. 

Theorem 5: Given GBGOP (M,s ,A,C, IC, C, Q in , 
Qout), if the cost function and all actions a € A can be 
polynomially computed, then determining if there is a solution 
SOL for the instance of the GBGOP s.t. for some real number 
k, \SOL\ <kis in-NP. 

Theorem 6: Given BMGOP (M,s , B, A, C, IC, k, c), if 
the cost function, benefit function, and all actions a £ A 
can be polynomially computed, then determining if there is 
a solution SOL for the instance of the BMGOP s.t. for some 

real number val, EA^appHsoL^) b i > val is in - NR 

As stated earlier, a GBGOP may also be viewed as an ab- 
ductive inference problem. Even though finding a solution (not 
necessarily optimal) to a GBGOP can trivially be conducted 
in PTIMBj, counting the number of solutions is #P-complete. 
This counting problem is difficult to approximate. 

Theorem 7: Counting the number of solutions to a GBGOP 
(under the assumptions of Theorem |5]l is #P-complete. 
Proof Sketch. The MONSAT problem [3| takes a set C of m 
clauses of K disjunct ed literals (no negation) over set L of 
atoms (size n) and outputs "yes" iff there is a subset of L 
that satisfies all clauses in C. We encode this into finding a 
GBGOP as follows: Q = {gx, . . . ,g m } - each predicate in Q 

4 Retum the set {(a^pi) e A X _M|a 4 (p i; ) n Bout = 9} 



corresponds to an clause in C (predicate gj corresponds with 
clause 4>j), Ai consists of a single point, p, A — {ai, . . . , a„} 
- each action in A corresponds to an element in L (action Oj 
corresponds with literal £j). Each Oj is defined as follows: 
a i(p) = {9j(p)\{^i} 1= $j}> C returns 1 for all action-point 
pairs, s = 0, IC = 0, C = n, <d in = \J gie g{9i(p)}, 
@owt = 0- Based on this PTIME reduction we show a 1-1 
correspondence to MONSAT. Hence, we can parsimoniously 
reduce the counting version of MONSAT (number of 
solutions) to the counting version of GBGOP (number of 
solutions). As the counting version of MONSAT is #P-hard 
by 0, we have shown that =^P-hardness of the counting 
version of GBGOP. As there is an obvious bound on the 
number of solutions to a GBGOP, and as the solutions are 
verifiable in PTIME, membership in jfP follows. □ 

Theorem 8: For e > 0, approximating the number of 
solutions to a GBGOP within a factor of 2 | - Ax - M|1 " is NP- 
hard. 

Follows from Theorem [7] and Theorem 3.2 of [0. □ 

Due to this issue with achieving a good approximation of 
the counting version, in this paper we shall focus only on 
determining a single optimal solution to a GBGOP - rather 
than all solutions. 

IV. Integer Programs for Solving GOPs 

In this section, we present an integer programming (IP) 
algorithms for both GBGOP and BMGOP which provide exact 
solutions. Given a GBGOP, the IP associates an integer-valued 
variable Xi with each action-point pair (ai,pi) G A x A4 
where di(pi) n & ut — 0- Intuitively, X; = 1 denotes that 
action a, is performed at point pi. 

Definition 4.1 (GBGOP-IP): Let set R = {{a u pi) € A X 
M\ai(pi)(l<d ou t = 0}. For each action-point pair (ai,pi) G R, 
create variable Xi £ {0, 1}. 

I «l 

min^X, (1) 

i = i 

s.t. X j> 1 Vj4 * G - s o (2) 

a j(Pj)\ A i£ a j(Pj) 

Y, c « • ^ < c (3) 

Xi<i V($^x)e^ so (4) 

The objective function minimizes the total number of action- 
point pairs. Constraint (01 ensures that every ground atom in 
Qi n (that does not appear in the initial state) is caused by 
at least one of the selected action-point pairs. Constraint (O 
enforces the constraint on cost. Constraint (@) ensures that the 
integrity constraints are satisfied. Next we present our integer 
constraints for a BMGOP where the IP associates an integer- 
valued variable X{ with each action-point pair (a,i,Pi) G A x 
M., and an integer-valued variable Yj with each ground atom 
Aj G Be — so. The intuition for the Xi variables is the same 
as in GBGOP-IP. 



Definition 4.2 {BMGOP -IP): For each action-point pair 
(a,i,pi) G A x A4, create variable Xj G {0,1}. For each 
A-i G Be — So create variable Yi G {0, 1}. 

IBjcI-M 

max Y b ' + Y bi ' Yl (5) 
s.t. J2 X i ^ Fl yA > £Bc~s (6) 

C r ■ X i < C (8) 

(a i ,p i )eAx M 

Y Xi < 1 V($ ^ x) G IC So (9) 

In the above IP, the objective function looks at each ground 
atom and sums the associated benefit if the associated Yi 
variable is 1 - meaning that atom Ai is true after the actions 
are applied. Constraint © effectively sets a Yi variable to 1 
if an action that causes Ai to be true occurs. Constraint (|7) 
enforces the cardinality requirement. Constraints |8]0 mirror 
constraints 04] of GBGOP-IP. The result below shows that 
a solution a to the above IPfl when restricted to the Xi 
variables, provides an immediate solution to the GOP. 

Prop. 4.1: Suppose T is a GBGOP (resp. BMGOP) and 
IP(T) is its corresponding integer program (GBGOP-IP, resp. 
BMGOP-IP). Then: 

1) If SOL is a solution to T, then there is a solution a of 
IP(T) such that a 2 {X = 1 | {a i:Pl ) G SOL}. 

2) If a is a solution to IP(T), then there is a solution SOL 
to T such that {X t = 1 | ((H,pi) G SOL} C a. 

As integer programming is NP-complete, any algorithm to 
solve a GOP using GBGOP-IP or BMGOP-IP using an IP 
solver will take exponential time. We note that for GBGOP- 
IP, the number of variables is fairly large - 0(|{(a,-,f>i) G 
A x M\cti(pi) n Qout = 0}|) variables and 0(\Q in - s \ + 
\IC So \ + 1) constraints. BMGOP-IP has even more variables - 
(though not exponential) - 0(|.M| • (|^4| + \G\)) variables and 
0(\M\ ■ \G\ + |7C So | + 2) constraints. However, BMGOP-IP 
has only packing constraints]! We also note the GBGOP-IP 
has both covering (>) and packing (<) constraints - another 
source of complexity. 

V. Correct Variable Reduction for GBGOP-IP 

The set of integer constraints for GBGOP has 0(|i?|) 
variables where R C A x Ai. We show how to correctly 
reduce the number of variables by considering only a subset 
of R - thereby providing a smaller integer program. Our 
intuition is that an optimal solution SOL is an irredundant 
cover of <di n meaning there is no subset SOL' C SOL that 
is also a solution. Hence, we can discard certain elements 

5 A solution to GBGOP-IP or BMGOP-IP is an assignment of values 
to variables that optimizes the objective function. Thus, a solution can be 
described as a set of equations assigning values to the variables Xi , Y 3 ■ . 

6 It is trivial to eliminate constraint |5]and re- write |5]as a non-linear objective 
function. 



of R that cannot possibly be in an optimal solution. First, 
for a given GBGOP T = (M,s Q , A, C, IC, C, 8 m , 6 „ t ), we 
introduce Q r (ap) = {$|(4> ^ X ) e IC So A (a,p) e $} 
and the set of ground atoms each action-point pair affects 
Aff[ n>p) = a 4 (pi) n (6j„ - (Q in n s ))- We can now define a 
reduced action-point set. 

Definition 5.1 (Reduced Action-Point Set): Given GBGOP 
T = (M, s , A, C, IC, C, 6 in , 9 out ) and set i? = {(a^p*) £ 
.4 x M|oi(pi) n Q out = 0}, we define reduced action-point 
set R* = {(ai,pi) £ i?| fl(a,j,pj) £ R s.t. 
( Cj < a) A (Ql jipj) C Qf Q<ij)i) ) A (Afff ai;Pi) C Afff„ JiPj) )} 

Example 5.1: Consider the campaign scenario 
last discussed in Example 12.51 Suppose the 
candidate wants to optimize the following GBGOP: 
r — ( ka a a c^ Sc ps") jr< a a c P9 n aw 

i \ J vl cpgn : cpgn 7 ^cpgn: w cpgn ' ±KJ cpgn: ^7 ^ in ' *V 

where each A £ Q c £ l 9n has the form exposure(p) where p 
is a point in one of the two dashed rectangles in Figure [TJ 
Note that as map M C pgn contains 187 points, _4| = 3, and 
Qout — 0, the cardinality of R is 561. By contrast, the set 
R* consists of only 7 elements, 1.2% of the size of R. Here 
R* = {(nor, (5, 4)), (nor, (5, 3)), (nor, (5, 2)), (nor, (10, 8)), 
(nor, (10, 7)), (nor, (10.6)), (appeah, (4, 3))} 

Intuitively, all elements in R* are preferable for membership 
in an optimal solution over R — R* as they cost less, result in 
the same changes to the state, and occur in the same or fewer 
integrity constraints. Set R* can be found in quadratic time 
with a naive algorithm - an operation that is likely dominated 
by solving or approximating GBGOP-IP. The next lemma says 
that R* must contain an optimal solution. any optimal solution 
to a GBGOP. This can then be used to correctly reduce the 
number of variables in GBGOP-IP. 

Lemma 5.1: Given GBGOP T = 

(M,so,A,C,IC,C,Qin,Q out ), for any optimal solution 
SOL C R, there is an optimal solution SOL' C 

Prop. 5.1: Suppose T is a GBGOP and LP(T) is its cor- 
responding integer program. We can create such a program 
with a variable for every element of R* (instead of R) and 
Proposition 14. 1 1 still holds true. 

VI. The BMGOP-Compute Algorithm 

While BMGOP-IP can solve a BMGOP exactly, doing so 
is computationally intractable. We now present an approxi- 
mation algorithm that runs in PTIME but provides a lower 
approximation ratio than proved in Theorem |4] First, we 
show that a BMGOP reduces to an instance of submodular 
maximization problem^ with respect to packing constraints. 
We then leverage some known methods |4) to solve such 
problems and develop a fast, deterministic algorithm to ap- 
proximate BMGOP with an approximation bounds. Given 
BMGOP T = (M,s ,B,A,C,IC,k,c), consider the objec- 
tive function in BMGOP-IP. We can write that function as 

1 Proof Sketch. We show this by proving that for any set W = SOLD(R — 
R*), there is some set W C R* - (R* n SOL) s.t. (SOL -W)UW' is 
also a solution. 

8 Suppose Z is a set. A function / : 2 Z — ¥ R is said to be submodular 
iff for all Z\, Z-z such that Z\ C Z2 and all z f Z2, it is the case that 
/(Zi U {z}) - /(Zi) > f(Z 2 U {z}) - f(Z 2 ), i.e. the incremental value 
of adding z to the smaller set Z\ exceeds the incremental value of adding it 
to the larger set Z2. Here, 1Z denotes the reals. 



a mapping from action-point pairs to reals. We denote this 
function (specific for BMGOP T) as f r : 2 AxM K+, where 
MS) = XUeappi^so) which has certain properties. 

fr(S)= b * ^ 

A 4 6appl(S,so) 

We now show that this function fr is submodular and has 
some other nice properties as well. 

Prop. 6.1: For BMGOP T, function f r is: (i) submodular, 
(ii) monotonic, i.e. Z\ C Z 2 -> fv{Z\) < /r(^) and (Hi) 
under the condition \/A { £ B c , h = 0, we have / r (0) = 00 
Proof Sketch. Consider S C S' C A x M and 
(a,p) i S'. We must show f r (S U {(a,p)}) - f r (S) > 
fr(S' U {(a,p)}) - fr(S')- Suppose, BWOC f r (S U 
{(o,p)}) - fr(S) < MS' U {(o,p)}) - /r(5')- Then, by 
Equation [TO] we have EA,ea W !(su{(a, P )}, s „)-a W i(s, ) 6 * < 

EA I 6a W i(S'U{(a, P )}, S „)-ap P i(S', S0 ) 6 - H ° WeVer ' b Y the 

definition of appl, we have appl(S U {(a,p)},so) — 
appl(S, so) 2 appl{S' U {(a,p)}, so) — appl(S' , s ), which is 
a contradiction. □ 

As our objective function is submodular, and constraints [TJ 
[9] are linear packing constraints, any instance of a BMGOP 
can be viewed as maximization of a submodular function wit 
linear packing constraints and hence, methods to solve such 
problems can be used here. The BMGOP-Compute algorithm 
leverages this idea and illustrated in Example 16.11 



BMGOP-Compute 

INPUT: BMGOP (M,s , B, A, C, IC, k, c) 
OUTPUT: SOL CAxM 

1) Set SOL — 0, 8 to be an infinitesimal, 

and set A = e 2 ' s ■ (2 + \IC SQ |). 

2) Set w' = 1/k and tu" = 1/c. For each ($» <-> x0 € -fC S0 , 
set id, = 1/(2 - 5). 

3) While fe-iu / +C-iu" + (2-a)-X) i Wi < A and SOL j^AxM 

a) Let (aj,pj) £ A x A4 — SOL have minimal 

w'+w"-c j +J2n(a :j , Pi ) ( =^ i m 

S^AieappKSOLUKaj , Pj )},s ) b i I E A; £ appl ( S O L . sq ) b i ) 

b) SOL = 50Lu{(a ilP3 )} 

c) Set w' = w' ■ X 1/k , w" = w" ■ A C J /C and for each 
integrity constraint i s.t. (a,j,pj) £ $i, set 

4) If SOL is not a valid solution then 

a ) lf T.A i eappl(SOL-{(a 3 ,p J )},s ) b ' ^ 

then SOL = SOL - {(a^,ft)} 

b) Else SOL = {(aj,Pj)} 

5) Return SOL 



Example 6.1: Following Example 12.51 Suppose 
the candidate wants to optimize BMGOP: 

o R A p( s cp9") Tfl q o\ T n 

cpgn: cpgn : L-tcpgn : -ST-cpgn: ^ cpgn 1 ±y - J cpgn: ^ )■ AI1 

this case, we will set 6 = 0.001. He wishes to find a set 
of 3 action-point pairs to optimize his exposure. BMGOP- 
Compute sets A = 22.14, w' = 0.33, w" = 0.50, and 

'Henceforth, we will assume this condition to be true. 



(M 



Wi — 0.50 in lines 1 and 2. In the first iteration of the 
loop at line [3] it finds the action-point pair that minimizes 
the quantity at line [3] is (appeah, (4, 3)) - which has the 
associated value 0.073. Note, other action-point pairs with low 
values are (appeal, (10, 7)) with 0.083 and (nor, (15, 6)) 
also with 0.083. It then adds (appeah, (4, 3)) to SOL and 
updates w' = 0.93, w" = 1.09, and wi = 2.35. On the next 
iteration, the BMGOP-Compute picks (nor, (15, 6)), which 
now has a value of 0.164. During this iteration, the value of 
(appeah, (10, 7)) has increased substantially - to 0.294, so 
it is not selected. At the end of the iteration, w 1 is updated 
to 2.611 and w" is updated to 2.364. As (nor, (15,6)) does 
not impact the lone integrity constraint, the value w\ remains 
at 2.354. In the third iteration, BMGOP-Compute selects 
(nor, (15, 9)) which has a value of 0.421. Again, the value of 
(appeal, (10, 7)) has increased - but this time only to 0.472. 
BMGOP-Compute re-calculates w' = 7.331, w" = 5.128 
and wi remains at 2.354. On the last iteration, BMGOP- 
Compute picks (appeah, (10,7)) as it has the lowest value 

- 0.942. After this fourth iteration, it updates w' = 20.589, 
w" = 11.124, and w 1 = 11.0861 - which now total to 42.799 

- exceeding A (22.14) - causing BMGOP-Compute to exit 
the outer loop. Now SOL has 4 elements, exceeding the 
cardinality constraint (as well as the integrity constraint). The 
checks done in line [4] remove (appeah, (10, 7)) from SOL 

- making the result feasible. BMGOP-Compute returns 
{(appeah, (4, 3)), (nor, (15, 6) j (nor, (15, 9))} which causes 
the benefit to be 45. 

Prop. 6.2: Suppose T is a BMGOP and SOL is the set 
returned by BMGOP-Compute. Then SOL is a solution to 

rE 

Next, we sho BMGOP-Compute runs in PTIME. 
Prop. 6.3: BMGOP-Compute runs in 0(k ■ \M\ ■ \A\ ■ 
\IC S „\) time. 

Proof Sketch. Clearly, the outer loop can iterate no more than 
k times. The inner loop iterates for each element of A x M - 
hence requiring time 0(|.M| • \A\). There are some additional 
operations that require 0(\IC Sn \) time, however, they are 
dominated under the assumption that \M.\ ■ \A\ >> |/C So |, 
which we expect in our application. □ 

The following important theorem states that BMGOP- 
Compute provides an approximation guarantee. Because of 
Theorem |4] and as BMGOP-Compute is polynomial, we 
know that this approximation guarantee cannot be as good 
as ^p- + e. The result leverages Theorem 1 . 1 of [4 1 together 
with the above theorems. By this result, the approximation 
factor of BMGOP-Compute depends on | IC So |. We illustrate 
this relationship, in Figure [2] For our target applications, we 
envision \IC So \ < 20. 

Theorem 9: Under the assumption that k, C > 2 — 6, 
BMGOP-Compute provides a solution within a factor of 
( 9+ | fC 1 |)i/(2-^) (where 6 is an infinitesimal) of optimal. 
Proof Sketch. BMGOP-Compute follows from Algorithm 1 

10 Here, SOL is not necessarily an optimal solution. 
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Fig. 2. |^C SQ | vs. approximation ratio. 

of H which optimizes a submodular function subject to m 
packing constraints within m i/ w where W is the minimum 
width of the packing constraints - defined as the minimum of 
the size of the constraint divided by the cost of an element. 
For constraint Q the W = k. For constraint [8] the W > C. 
We can replace constraint [9] with: 2~Z(a - — 2 — <5 

V(<i>j <r^ Xj) S IC So which maintains correctness as two 
variables to set to 1 and exceeds 2 — 6. The new constraint 
has width 2 — 5, which, is the minimum. We then apply 
Theorem 1.1 of ||4). □ 

Discussion. We note that while a BMGOP reduces to the max- 
imization of a submodular or linear function wrt linear packing 
constraints, there are other algorithms available besides the 
multiplicative update algorithm of [4|. However, we feel that 
this is likely the best approach for several reasons that we list 
below. 

1) The approximation ratio achieved by the multiplicative- 
update algorithm matches the best approximation ratio 
achievable for maximizing a linear function wrt linear 
packing constraints (see 0), hence, it is unlikely that a 
better approximation ratio can be achieved using such a 
technique. 

2) Other methods (such as those presented in 0) require 
solving a relaxation of the associated MILP. In our case, 
such an operation would take 0((\M\ ■ (\A\ + |^|)) 3 ' 5 ) 
time (as a consequence of the number of variables in 
BMGOP-MILP and the results of 0). This is signif- 
icantly more expensive than the 0(k ■ \A4\ ■ \A\) of 
BMGOP-MU (see Proposition l63l If the map, M is 
very large, solving a relaxation of BMGOP-MILP may 
be unrealistic on most hardware. 

3) The algorithm BMGOP-MU is totally deterministic, 
which allow us to avoid the issue de-randomization. 

4) The algorithm BMGOP-MU is guaranteed to provide a 
solution that meets constraints |7]|9] - as opposed to only 
meeting them probabilistically. 

VII. Related Work and Conclusions 

Though spatial reasoning has been studied extensively in AI 
0, 10, 13, iflOl . many of the paradigms that have emerged 
for such reasoning are qualitative in nature. Such qualitative 
spatial reasoning efforts include the influential region con- 
nection calculus for qualitative reasoning about space. There 
has also been work on quantitative methods for reasoning 
about space ifTTI which contains articles on spatial reasoning 
in the presence of uncertainty using both logical and fuzzy 



methods. Spatial reasoning with quantitative information has 
been studied extensively in image processing fl2l . fl3l . 

However, unlike this vast body of work, this paper focuses 
on a different problem. Suppose we are dealing with a map 
A4, a cost function C, a set A of possible actions, a bound 
on the cost C, and a bound on the number of actions we can 
take, what set of actions should be taken so as to optimize a 
given objective function. Two versions of this problem are 
studied in this paper - GBGOP and BMGOP which differ 
in what they optimize. Both problems are proved to be NP- 
hard (NP-complete under realistic assumptions) and we further 
prove that the number of solutions to GBGOP is #P-complete. 
We also find limits on approximating an optimal solution to 
BMGOP and GBGOP (in PTIME) under accepted theoretical 
assumptions. We develop integer programming formulations 
of both problems and then present a way of simplifying the 
IP for GBGOP. We further present the BMGOP-Compute 
algorithm for BMGOP and show that it is polynomial and has 
a guaranteed approximation ratio (though not high enough to 
contract the NP-hardness result). 

VIII. Conclusion 

In this paper, we introduced "geopspatial optimization prob- 
lems" or GOPs that aide the user in taking certain actions over 
a geographic region. We showed these problems to be NP-hard 
and provided integer constraints. For the goal-based variant, 
we correctly reduce the number of variables. For the benefit- 
maximizing variant, we provide an approximation algorithm. 
In future work, we look to implement this framework and 
explore methods to achieve further scalability, as well as uti- 
lize geo-located social network data to establish relationships 
among locations in order to better implement action-point pairs 
and integrity constraints. 
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